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ABSTRACT 


We study the rotational distortions of the vacuum dipole magnetic field in the 
context of geometrical models of the radio emission from pulsars. We find that 
at low altitudes the rotation deflects the local direction of the magnetic field by 
at most an angle of the order of rj;, where r n = r/Ri c , r is the radial distance and 
R lc is the light cylinder radius. To the lowest (ie. second) order in r n , this distor- 
tion is symmetrical with respect to the plane containing the dipole axis and the 
rotation axis ((fi,/2) plane). The lowest order distortion which is asymmetrical 
with respect to the (Q, p) plane is third order in r n . These results confirm the 
common assumption that the rotational sweepback has negligible effect on the 
position angle (PA) curve. We show, however, that the influence of the sweep- 
back on the outer boundary of the open field line region (open volume) is a much 
larger effect, of the order of r]J 2 . The open volume is shifted backwards with 
respect to the rotation direction by an single S ov ~ 0.2 sin ar ]/ 2 where a is the 
dipole inclination with respect to the rotation axis. The associated phase shift 
of the pulse profile A<p ov ~ 0.2 tv ' 2 can easily exceed the shift due to combined 
effects of aberration and propagation time delays (« 2 r n ). This strongly affects 
the misalignment of the center of the PA curve and the center of the pulse profile, 
thereby modifying the delay-radius relation. Contrary to intuition, the effect of 
sweepback dominates over other effects when emission occurs at low altitudes. 
For r n < 3 • 10 -3 the shift becomes negative, ie. the center of the position angle 
curve precedes the profile center. With the sweepback effect included, the mod- 
ified delay-radius relation predicts larger emission radii and is in much better 
agreement with the other methods of determining r n . 


Subject headings: pulsars: general — polarization — radiation mechanisms: non- 
thermal 
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1. Introduction 


There are independent observational arguments which imply that the pulsar radio emis- 
sion occurs in a form of a narrow beam centered (or roughly centered) on the magnetic 
dipole axis (eg. Radhakrishnan & Cooke 1969; Lyne & Manchester 1988; Rankin 1990; 
Rankin 1993). It is commonly believed that the emission region associated with the beam 
does not extend beyond the region of open field lines (hereafter called “open volume” ) which 
cross the light cylinder of radius Ri c = c/Q (c is the speed of light and Q. = Q,z is the angular 
velocity of pulsar rotation). The angular size of the open volume at (small) radial distance 
r is equal to 9 r ov ~ {r / R \ c ) 1//2 and the cone formed by tangents to magnetic field lines at the 
rim of the open volume has angular radius of ^ 1.5#£ v . The radial distance of the emission 
region has not been established so far: both a high-altitude emission region extending over 
a small fraction of 0£ v , as well as a low-altitude emission region which fills in a much larger 
fraction of 9 r ov may be responsible for the same shape of the radio beam. 

The sweepback effect was first investigated in detail by Shitov (1983) who considered it 
to explain the observed dependence of radio luminosity of pulsars as a function of period. He 
estimated the magnitude of the rotational distortions of the magnetic field from the torque 
responsible for the observed slowing down of pulsars. He found that at moderate altitudes 
within the open volume, “near” the dipole axis, the direction of the distorted magnetic field 
deflects from the direction of the pure (ie. static shape) dipole barely by an angle 
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(1) 


where a is the dipole inclination with respect to the rotation axis. Since that time the 
sweepback effect has been recalled to explain asymmetries observed in some pulse profiles 
(eg. Gil 1983). 

In 1985 Shitov incorporated the sweepback effect into the model of pulsar position angle 
curves proposed by Radhakrishnan & Cooke (1969) and showed that the sweepback results 
in a lag of the profile center (measured as the midpoint between the outer edges of the 
pulse profile) with respect to the center, or the “inflection point” of the position angle curve. 
Shitov emphasized that the lag of the profile center was a sum of two effects: not only the 
center of the PA curve is shifted toward earlier phases (with respect to the nondistorted 
case) according to the eq. (1), but also the center of the open volume is displaced backwards, 
which contributes to the total effect. 


In most of subsequent investigations, however, the sweepback has been neglected, mainly 
on the basis of eq. (1). Blaskiewicz et al. (1991, hereafter BCW91) proposed a relativistic 
model of pulsar polarization which took into account two important effects overlooked by 
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Shitov: the presence of the corotational acceleration and the aberration effect. An excellent 
result of their work was the “delay-radius” relation, according to which the center of the PA 
curve lags the profile center by 

bcw ~ 4-r— - rad, (2) 

ttlc 

where r is the radial distance of the radio emission. With no dependence on viewing geometry 
parameters (like the dipole inclination a, or the viewing angle Cobs between the rotation 
axis and the observer’s line of sight), their relation appears to provide a powerful method of 
determining r. Equally important, the delay-radius relation depends neither on the observed 
width of the pulse profile W nor on the separation between the conai components in the 
pulse profile. Therefore, the altitudes of radio emission provided by eq. (2) may serve to 
determine which magnetic field lines are associated with the outer edge of the profiles and 
which field lines correspond to the maxima of conai components (Mitra & Rankin 2002; Dyks 
et al. 2004a). Von Hoensbroech & Xilouris (1997) used the delay-radius relation to probe 
the radius-to-frequency mapping at high radio frequencies. 

Gangadhara & Gupta (2001) proposed another relativistic method of estimating radio 
emission altitudes for pulsars with both core and conai components. By considering the 
effects of the aberration and the propagation time delays they showed that the core com- 
ponent lags in phase the midpoint between the maxima of conai components, if the core 
originates from lower altitudes than the cones, and if the cones are axially symmetric around 
the core in the reference frame corotating with the star (CF). Dyks et al. (2004a) revised 
their method and showed that the phase shift between the core component and the pairs of 
conai components is equal to 

^■^drh ^ 2-rr— rad (3) 

-Klc 

which provides another method for determining r without information about viewing geome- 
try (nor W). As in the case of the delay-radius relation,- the above formula holds only for the 
magnetic field which is symmetrical about the (D, p) plane (where p is the magnetic moment 
of the pulsar magnetic field), at least as long as one associates the assumed symmetry of the 
core-cone system with the geometry of the underlying magnetic field. 

Given that the methods provided by eqs. (2) and (3) are based on a measurement of tiny 
phase shifts (of magnitude usually being a small fraction of one degree) they are extremely 
sensitive to the assumed geometry of the magnetic field. The latter was taken to be a 
dipole of static shape, with no rotational distortions. It is therefore important to study the 
influence of the sweepback on these methods. Another argument for the study is provided 
by the unacceptably low values of r which are often being derived with the BCW91 method: 
as found in BCW91, the “delay radii” rjei implied by their method (eq. 2) are often smaller 
(in some cases by an order of magnitude — see fig. 29 in BCW91) than the geometrical radii 
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r geo determined with the traditional geometrical method based on the measurement of profile 
widths (Cordes 1978; Gil &; Kijak 1993; Kijak & Gil 2003). This poses a real problem for the 
BCW91 method, because the geometrical radii, in the absence of strong refraction effects 
(Lyubarski & Petrova 1998), should be considered as lower limits of r (Dyks et al. 2004a). 
The reason for this is that we most probably miss the emission from the outer parts of 
the open volume, either because this region is inactive, or because of the limited sensitivity 
threshold (or because of both). Although strong irregularities of the radio beam could result 
in r<iei < r geo in some cases, on average the opposite trend should be observed, but it is not 
(eg. Mitra & Li 2004). Given that BCW91 “cut away” the outer wings of pulse profiles by 
assuming the lowest intensity points at the high 10% level, the “center of mass” of points 
in their Fig. 29 should be located clearly above the diagonal of perfect agreement. In this 
paper we argue that the rotational distortions of the static shape dipole are responsible for 
a large part of this problem. Among the other effects which work in the same direction, one 
we consider is the variation of emission altitude across the pulse profile, mentioned already 
by BCW91 and described in detail in Dyks et al. (2004a). 

Recently Kapoor & Shukre (2002) considered the aberration effect and the rotational 
sweepback to investigate the relative locations of core and cone components in the pulsar 
magnetosphere. Although included in the model, the sweepback is again estimated with the 
help of eq. (1). Being aw r are of the limitations of Shitov’s estimate, the authors emphasized 
the need for derivation of a more advanced formula describing the rotational distortions of 
the magnetosphere. They noted that a proper derivation “should make use of at least the 
magnetic field given by the full Deutsch solution (Deutsch 1955)” . 

Such an estimate based on the Deutsch solution was done by Arendt & Eilek (1998), 
who concluded that the rotation distorts the magnetic field by a magnitude of the order of 
r/R\ c . Being much larger than the Shitov’s estimate, this distortion would strongly affect 
results in BCW91, GG2001, Hibschman & Arons (2001, hereafter HA2001), and Dyks et 
al. (2004a). On the contrary, HA2001 noted that the leading terms in the difference between 
the Deutsch field and the rigidly rotating static-shape dipole are of the order of ( r/R \ c ) 2 . 
Recently, Mitra &. Li (2004) emphasized that on the theoretical side there is a great need to 
develop and understand the details of the sweepback effect. 

In this paper we investigate the rotational distortions of the pulsar magnetic field assum- 
ing the approximation of the vacuum magnetosphere. The twofold nature of the sweepback, 
first noticed by Shitov (1983) will be highlighted, and limitations in applicability of eq. (1) 
will be clarified (Section 2). The significance of the sweepback for the relativistic model of 
pulsar polarization will appear to be much larger than previously thought, which will have 
serious consequencies for the delay-radius relation (eg. modification of eq. 2, Section 3). 
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2. The rotational distortions of the dipolar magnetic field 


We follow previous investigators (Deutsch 1955; Shitov 1983; Barnard 1986; Romani & 
Yadigaroglu 1995, hereafter RY95; Cheng et al. 2000, hereafter CRZ2000) in assuming that 
the magnetic field surrounding the neutron star (NS) may be approximated by the vacuum 
rotating dipole. As in Barnard (1986), RY95, and CRZ2000, we assume that outside the 
NS the field is the same as of the star-centered point dipole, ie. we neglect the near-surface 
modifications of the magnetic field by the conducting sphere of the neutron star, derived 
by Deutsch (1955) (see Yadigaroglu 1997). Hereafter, the magnetic field will be called a 
“retarded dipole” and will be denoted by B Kt . In Appendix A we give the cartesian and the 
spherical components of B ret (eqs. Al - A3 and A15 - A17, respectively). 

We want to estimate how much the rotational sweepback distorts the magnetic field at 
low altitudes (r <C R\ c ). One measure of this is the difference between the retarded magnetic 
field B Tet and the magnetic field of the static-shape dipole B^. The components of B st can 
be calculated with the help of eqs. (Al - A3) taken in the limit of r n <C 1 (ie. with the ratio 
r n set equal to zero). We define the difference as: 

AB — S ret — B s t- (4) 

In all formulae we assume that both the retarded dipole and the static-shape dipole are 
associated with the same magnetic moment j I, which at the time t = 0 is in the (x, z) plane 
(time t is measured in the Lorentz frame in which the neutron star’s center of mass is at 
rest). Thus, at any instant /7 st = /2 ret = /Z, where 


jl = p.(sin a cos f It x+ sin a sin f U y + cos a z ). 


(5) 
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At any point which corotates with the magnetosphere the components of B ret in the inertial 
observer frame (IOF) do not depend on t (see eq. A4). Therefore, one is allowed to choose 
any convenient value of i in the IOF. We take i = 0 (/Z in the (x, z) plane) and constrain our 
discussion to the half of the magnetosphere with positive values of x. The positive values of 
the y coordinate then correspond to the leading part of the magnetosphere and the negative 
y correspond to the trailing part. In cartesian coordinates and for t = 0 the difference is: 


A B x 
A B y 
A B z 


V . 

— sin a 




sin a 




sma 


\ (* 2 + r2 ) r l + °( r l) 

\ *y r l - \ x 2r l + 0{r 4 n ) 

\ xz r 2 n + 0(r*) . 


( 6 ) 

(7) 


(8) 
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In agreement with the remark of Hibschman &; Arons (2001), the leading terms are second 
order in r n . The second order terms of A B x and of A B z do not depend on y. The second 
order term of APj, is odd function of y (A B y (y) — —A B y {-y) + 0(r®)). These features 
are important because the symmetry of any vector field B = (B x , B y , B z ) with respect to 
the (Q, jT) plane requires the following relations to be satisfied: B x (x,y,z ) = B x (x,—y,z), 
B y (x,y,z) = -B y (x,-y,z), and B z (x,y,z ) = B z (x,—y,z) (ie. the B x and B z components 
must be even, and the B y component must be odd in y). 

The angle k between P ret and P st , to the order r^, is given by: 

kz* (hrt + hr 5 n + f 3 rt) 1/2 <rl, (9) 


where the functions fi, / 2 , and / 3 , given in Appendix A, depend on x, y, z, and on the 
inclination angle a (but not on r n ) and in general have magnitude of the order of 1, except 
from special locations in the magnetosphere which we discuss below. 


Thus, the rotation causes the magnetic field to deviate from P st by an angle which 
at most is second order in r n . Along the magnetic dipole axis 2 of an orthogonal rotator, 
however, (ie. for a = 90° and (x,y,z) = (r, 0,0)), one obtains /i = 0, / 2 = 0 and / 3 = 9 -1 
so that k is third order in r n : 



( 10 ) 


in partial agreement with the estimate of Shitov (1983). Beyond the orthogonal dipole axis, 
however, as well as at the dipole axis of non-orthogonal rotator, fi^O and k may be much 
larger. On the (Q,jT) plane /2 = 0 (because f 2 oc y, see eq. A-5). Beyond the (Q,j2) plane 
(/i 7^ 0, / 2 0), the first two terms in eq. (9) dominate and give: 


K 




h r 3 

1/2 n* 


2/; 


(ii) 


Let us estimate the 'angles K\ and « t for two points P\(x\, y\, z{) and P t (x t , y t , z t ) located 
symmetrically on both sides of the (Q, ft) plane. Let us consider the particular case in which 
the points lie in the plane of rotational equator, close to the rim of the open volume of 
orthogonal rotator, ie. x\ = x t = r( 1 — r n Y / 2 , y\ = —y t = rr^ 2 , and z\ = z t = 0, with the 
positive value of the y coordinate corresponding to the point Pi on the leading side of the 
open volume, and the negative y for the trailing point P t (see Fig. 1). Then eq. (A5) gives 
/i = 4 -1 r n , / 2 = ±3 _1 ry 2 , which results in k\ » 2~ l rn 2 + 3 -1 rjj and K t w 2 _1 r^ 2 — 3 -1 r^. 
Both angles are considerably larger than the distortion of the magnetic axis given by eq. (10). 


2 By the “magnetic dipole axis” we understand the straight line containing the magnetic moment p. In 
the case of the retarded dipole, the axis cannot be associated with any magnetic field line. 
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Eqs. (6 - 8) imply, however, that up to the second order in r n , the rotational distortion 
of B st is symmetrical with respect to the (ft, jl) plane (because A B x is even, and the leading 
term of A B y is odd in y). This (approximate) symmetry, shown in Fig. 1, implies that 
beyond the dipole axis of the orthogonal rotator the angle k between B iet and B sl provides no 
information about the magnitude of the asymmetry of £ ret with respect to the (ft, jl) plane. 

To estimate the asymmetry for points located beyond the (ft, jl) plane, one must there- 
fore use the difference between azimuths of B vet at the points Pi and P t : 


A0i_ t = \4>t\ -4>\~ 


Bret ,y (Pt) 

B r et,x{P t) 


BretjA 1) _ sin a: 3 
iW^)~3r 3 B st>x D ’ 


( 12 ) 


where Q\ is the azimuth of P ret at the point P\{x\, yi, zi), 4>t is the azimuth of B ret at 
Pt{x h -y h z x ) (see Fig. 1), and 


Bgt,x — 


E 

3 


3 — cos a + 
r r 



sin a . 


(13) 


Equation (12) clearly demonstrates that Ap_ t is third order in r n , ie. the rotation induces 
the asymmetry of B ret with respect to the (ft, jl) plane with magnitude of the order of r 3 . 
Equation (12) is not useful in the immediate vicinity of the (ft, jl) plane (nor at the plane 
itself), because (j>\ changes sign to negative (ie. B ret is parallel to the (ft, jl) plane) for locations 
with the tiny azimuth 


4>\\ 



2 r s /££ JL + f^V 1 * £ 

9 n \ r r tan a r 2 ) 9 


(14) 


(ie. within a narrow region on the leading side of the (ft, jl) plane ^ret ,y < 0 and eq. (12) 
gives a sum rather than a difference of azimuths). To estimate the asymmetry on the (Q, jl) 
plane one can use the difference of azimuths of B st and B ret at a given (the same) point: 


Afis — r — 4*st 4 * ret 


P 

B t 


st ,y 


R 


'st, i 


B 


^ret,y 

ret,x 


U r l + ? 


jJi sin Gi 0 


3 r 3 B s t, x 


n’ 


(15) 


where the function f 4 (x, y, z, a) is given in Appendix A (eq. A14). Beyond the (ft,/?) plane 
the first term in this equation dominates (ie. A </> s _ r ~ f 4 rf) but it is symmetrical with 
respect to the (ft, jl) plane (ie. odd in y, / 4 oc y). Therefore, just like k given by eq. (9), 
A</> s _ r provides no estimate of the rotational asymmetry there. On the (ft, jl) plane / 4 = 0 
and A4> s -r does measure the asymmetry which is of the order of r 3 . Along the dipole axis 
P st ,i = 2 p sin ct/r 3 so that A 4> s -r — 3 -1 r 3 , independent of a. 

Eqs. (9), (12), and (15) can be summarized as follows: the rotation changes the compo- 
nents and the direction of the dipolar magnetic field by ~ r„. To the order of r% this change, 
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however, is symmetrical with respect to the (£2, jl) plane. The asymmetrical change of the 
magnetic field direction has much smaller magnitude of the order of rjj and is given by the 
second term of eq. (15). 

An immediate consequence of eq. (9) is that with accuracy of the order of r n , the rotation 
does not affect the shape of the position angle curve which depends on the local direction of 
J3 ret . This allows one to neglect the influence of the sweepback on the position angle curve 
as long as only the first order effects in r n are considered. Given the tiny magnitude of the 
rotational asymmetry as defined by eq. (12), investigators often neglect its influence on the 
shape of the pulse profile as well (eg. BCW91, GG2001, HA2001). As we show below, this 
is not justified. 


2.1. The rotational distortions of the open field line region 

In the method of BCW91, the center of the pulse profile is most efficiently measured 
as a midpoint between the outer edges of the profile. Therefore, the method is based on 
the assumption that the outer boundary of the open volume is symmetrical with respect to 
the (ft, jl) plane. Due to the complexity of the magnetic field lines in the retarded case, we 
determine the outer boundary of the open volume numerically, by finding the magnetic field 
lines which are tangent to the light cylinder and never cross it. The method is described in 
detail in Dyks et al. (2004b). 

Thick solid lines in Fig. 2 present the transverse shape of the open volume at low 
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altitudes (r n C 1), calculated for the retarded magnetic field B vet . More precisely, they 
represent the crossection of the outer boundary of the open volume with a sphere of radius 
r = 0.01-Ric centered at the neutron star. Different panels correspond to different dipole 
inclinations a. The magnetic moment jl in all panels emerges perpendicularly from the 
page at the point (x m ,y m ) = (0,0). The thin circles have radius equal to r0 r ov = rry/ 2 , 
and are centered at the (0,0) point to guide the eye in assessing the asymmetry of the 
open volume around jl. Given the small difference between the local direction of B ret and 
B st (eq. 9), one may regard each panel to be permeated by the magnetic field of the static 
dipole with the straight magnetic field line emerging from the (0, 0) point toward the reader. 
The field is symmetric relative to the (f2, /T) plane, ie. with respect to the vertical line of 
x m — 0 (to be imagined in each panel of Fig. 2). In the course of corotation, the contours 
of the open volume outer boundary move to the left in Fig. 2, ie. x m < 0 correspond to 
the leading, and > 0 to the trailing side. An observer’s line of sight cuts the contours 
horizontally, moving left to right. The dotted lines present the shape of the open volume 
for the static dipole (eg. Roberts & Sturrock 1972; Biggs 1990; Beskin, Gurevich & Istomin 
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1993; Shukre & Kapoor 1998). We emphasize that the definition of the open field line region 
based on the light cylinder radius is inherently inconsistent with the static dipole geometry. 
The ovalness generated by the inclination of the dipole is therefore a premature refinement, 
which practically disappears as soon as the rotation is taken into account. 

We emphasize that our Fig. 2, calculated for t — 0 in eqs. (A1 - A3), differs significantly 
from analogous figures shown in Arendt & Eilek (1998) and in CRZ2000. The location of the 
retarded polar caps in their figures corresponds to the magnetic moment ft rotated by the 
angle f IR^/c with respect to the (x, z) plane (ie. they assume t = R^/c in eqs. Al -A3). At 
the same time, however, they assume t = 0 (/! in (x, z) plane) to position the polar caps for 
the static case. Therefore, their figures do not inform us what is the relative position of the 
static and retarded caps in phase - only the caps’ shapes can be compared. (To enable this, 
their retarded caps would have to be derotated by £lR ns /c with respect to the static caps.) 
Also, note that the components of the magnetic field for the static-shape dipole given in 
CRZ2000 (eqs. A1-A3 therein) are for t = 0 whereas the components for the retarded dipole 
(eqs. B2-B4 in CRZ2000) are for t — R ris /R\ c . Their difference does not give eqs. (6 - 8): 
by overlooking this misalignment of the dipoles Arendt & Eilek (1998) incorrectly estimated 
that the rotational distortions of the magnetic field are of the order of r n . Another issue 
related to Fig. 2 is that for moderate dipole inclinations (a ~ 40° — 50°) the contours possess 
a notch rather than a discontinuous “glitch” suggested in Arendt & Eilek (1998) (see Sec. 2.1 
in Dyks et al. 2004b for details). 

The following important conclusions can be drawn from Fig. 2: 1) the open volume is 
strongly asymmetric with respect to the (fl, /Z) plane; 2) the magnitude of the asymmetry 
depends on y m and thus on the impact angle /3; 3) regardless of the value of (3 the pulse 
window associated with the outer boundary of the open volume is always shifted backwards, 
ie. toward later phases; 4) the magnitude of this rotational asymmetry is very large, much 
larger than the local changes of the direction of the magnetic field caused by the rotation 
(eq. 9), and even larger than r n (for large a the “retarded contours” are (on average) shifted 
toward later phases by ~ O.20” v = 0.2fn /2 = 0.02 to be compared with r n = 0.01). 

This numerical result implies that none of the previous estimates (neither eqs. 12, 15, 
nor the Shitov’s formula 1) provide a reliable measure of the rotational distortion of the 
open volume shape. The reason for this is that the boundary of the open volume is not only 
determined by the local (ie. low- altitude) direction of B ret , but also (and most importantly) 
by the geometry of B ret near the light cylinder, where r n ~ 1. At R\ c all “higher order” effects 
become comparable in magnitude to the lowest order effects, in the sense that r™ ~ 1 for any 
m. With the strength of the rotational distortions being very large at R\ c , a very different 
set of magnetic field lines is picked up as the last open field lines which form the boundary 
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of the open volume. This “retarded” boundary is highly asymmetrical with respect to the 
(0, jl) plane (Fig. 2). The magnitude of this asymmetry has little to do with the low-altitude 
rotational distortions as estimated with eqs. (12), (15) and with Shitov’s formula, because 
the low-altitude crossection of the open volume boundary is an image of the strong near-R\ c 
distortions projected through the continuity of the magnetic field lines. Hereafter, we will 
refer to the asymmetrical distortion of the open volume with the terms “backward shift” or 
“displacement” of open volume. 

Another numerical result is that the contours shown in Fig. 2 are (with high accuracy) 
the same for any radial distance r n , as long as r n <C 1, and as long as their size is normalized 
by r6 T ov = rr„ /2 , as in Fig. 2. Contours calculated eg. for r n = 10~ 3 or r n = 0.1 look exactly 
the same as those shown in Fig. 2. This means that for r n <l the overall (ie. averaged over 
the impact angle j3) angular displacement of the open volume S ov is a fixed fraction of the 
angular radius of the open volume 9 T 0V \ 

5 0V ~ ad r ow « ar ]/ 2 (16) 

with a ~ 0.2 sin a. The dependence a oc sin a has been determined by noting that 5 0v 
(measured from the star center) decreases with a (Fig. 2), whereas the corresponding phase 
shift A</> ov « 6 0v / sin a of the pulse profile does not (see Fig. 3, the discussion in Section 3.1, 
as well as eq. A8 in Dyks et al. 2004a). 

Equation (16) implies that the rotational displacement of the open volume has magni- 
tude comparable to the combined effects of aberration and propagation time delays (hereafter 
APT effects, of magnitude 2r n , see Dyks et al. 2004a) for r n ~ a 2 / 4 ~ 0.01, which is quite 
typical estimate of radio emission altitudes (eg. Gupta & Gangadhara 2003; Kijak &; Gil 
2003). For r n <§; 0.01 the backward shift of the open volume dominates over the APT 
effects. For r n 0.01 the APT effects dominate over the open volume shift. 


2.2. Twofold nature of the rotational distortions 

We find, therefore, that the nature of the rotational distortions at low altitude is twofold: 
in addition to the famous (but negligible) asymmetrical distortion of the local magnetic field 
direction (of magnitude ~ r 3 , eqs. 12 and 15) the rotation shifts the open volume backward 
by a much larger amount of ~ r\! 2 . The low-altitude shift of the open volume is not caused 
by the distortions of the shape of magnetic field lines at low altitudes - locally their shapes 
are pretty much the same as those of the static dipole. It is due to the strong near-f?i c 
distortions. 

Although Shitov (1985) in his analysis of the phase shift between the pulse profile 
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center and the center of the position angle curve neglected the APT effects, he did include 
the backward displacement of the open volume and did emphasize the twofold nature of the 
rotational distortions (see his fig. 2). However, he has not provided any simple estimate of the 


open volume displacement (like eq. 16). Therefore, in most studies following his work only 
the tiny local deflection of B has been considered (eqs. 12, 15 and Shitov’s estimate), usually 
only to infer that the rotational distortions are negligible in generation of any asymmetry in 
a pulse profile. 


3. Implications for the relativistic model of pulsar polarization 

Implications of the rotational displacement of the open volume for the relativistic model 
of pulsar polarization are profound, because this effect is lower order (~ rl /2 , eq. 16) than 
the effects considered so far (aberration ~ r n , propagation time delays ~ r n ). Although the 
PA curve is practically unaffected, the center of the pulse profile to which the PA refers is 
considerably displaced. 

Let us define a phase zero as a moment at which an observer detects a light signal emitted 
from the neutron star center when Jx was in the (£T, /2) plane. As discussed in Dyks et 
al. 2004a, the total phase shift between the center of the position angle curve and the 
pulse profile center can then be separated into two components: the shift of the center of the 
PA curve by Ad>p A « 2 r n towards later phases with respect to the zero phase, and the shift 
of the pulse profile center toward earlier phases by A<p pf with respect to the zero phase. Had 
the boundary of the open volume been symmetrical with respect to the (Q, /x) plane (as in 
the case of the static-shape dipole) the profile center would be shifted forward in phase by 
Ad>p f — 2r n which would result in the total shift of 4r n as initially predicted by BCW91. 
Due to the backward displacement of the open volume given by eq. (16), however, the forward 
profile shift is decreased by A (p ov = Fr]J 2 (with F ~ a / sin a) so that Ad> pf = -(2r n - FrlJ 2 ). 

Therefore, the delay-radius relation of BCW91 (eq. 2) becomes: 

A$ = A(p PA + (-Atfpf) « 2r n + (2r n - FFJ 2 ) (17) 

with F in general being a complicated function of a, £ and r n . We emphasize that this formula 
holds only when the detectable radio emission fills in the open volume. If the beam’s outer 
edge is irregular, the formula may give erroneous results for a particular object, however, it 
is still useful for statistical studies. The complicated form of F results from the complicated 
shape of the open volume boundary (Fig. 2) which implies nontrivial dependence of F on 
the impact angle and thereby on £. Since the sign of the impact angle 8 = ( — a provides no 
information about whether the viewing trajectory is poleward or equatorward (in the sense 
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defined in Everett & Weisberg 2001) hereafter we will use the angle 


f p, if a < 90° 
\ ~P, if a > 90° 


(18) 


which is negative/positive for poleward/equatorward viewing geometry regardless of whether 
a > 90° or not. Changes of F as a function of r = P x /@ b ^ Px/( l-Sr^ 2 ) are illustrated in 
Fig. 3 for the same angles a as those in Fig. 2. The functions F(r) were calculated for 
r n = 0.01, however, they change little with r n , as long as r n <C 1. Fig. 3 shows that F is 
confined to the rather limited range of 0.1 — 0.4 for any combinations of a and C- Also, F is 
always positive which implies that the displacement of the open volume results in a smaller 
phase shift A$ than predicted by the original delay-radius relation (eq. 2). The radio 
emission radii provided by the original delay-radius relation are therefore underestimated by 
a factor which may be very large for some parameters. 


3.1. The misalignment formula 

The delay-radius relation which includes the rotational distortions of the open volume 
(in the vacuum approximation) becomes: 

A«F«4r n -F(o,C,r n )rJ;/ 2 (19) 

and it is plotted in Fig. 4 for a few values of F equal to 0.1, 0.2, 0.3, and 0.4 (solid, dotted, 
dashed, and dot-dashed curve, respectively). (For convenience, A$ on the vertical axis is 
in degrees, whereas eq. 19 gives A<F in radians.) The thick dashed line presents the original 
delay-radius relation of eq. 2. For r n < F 2 / 16 (eg. for r n < 10 -2 for F = 0.4, dot-dashed 
line), ie. for small emission radii, the phase shift A$ becomes negative (ie. the center of the 
PA curve precedes in phase the profile center), which is a new feature in comparison with the 
original delay-radius relation which always predicted the delay of the PA curve with respect 
to the profile. Since eq. (19) predicts that the center of the PA curve may either precede 
or lag the center of the profile, it will be referred to as the “misalignment” formula. For 
large radii r n > 10 -2 the formula always predicts positive A$ and, for increasing r n , slowly 
converges to the original delay-radius relation. 

For any A$, the radio emission radii estimated with the original delay-radius relation 
underestimate those implied by eq. (19) (see Fig. 4). For A<E> ~ 1°, the delay-radius relation 
underestimates r given by eq. (19) by a factor of 1.5 — 4 (depending on F). For A4> ss 0.3°, 
0.1°, and 0.01° the underestimate factor is in the range of 2 — 9, 3 — 22, and 16 — 220, 
respectively. For A4> ~ 0 the underestimation factor is formally infinite. In the absence 
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of effects described in Section 4, the negative values of the shift AT cannot be lower than 
-F 2 / 16 > -0.57°. 

Unlike the original delay-radius relation, the misalignment formula predicts that the 
implied emission radii depend on the viewing geometry (especially on J3 X ) and that this 
dependence cannot be separated from the dependence on r n : even for fixed a and /? x , when 
the line of sight probes the magnetosphere deeper and deeper (ie. when r n decreases in 
Fig. 4), the absolute value of the parameter r = /? x /(1.57v ) increases — the line of sight 
cuts through the open volume more peripherally. This departure from r = 0 in Fig. 3 implies 
that the value of F changes (sometimes abruptly) with varying r n , ie. for fixed a and /? x the 
value of F is not fixed — it depends on r n . 

The complicated behaviour of the misalignment formula (19), is exemplified in Fig. 5, 
which presents AT as a function of r n calculated numerically for a = 45° and three values 
of Cobs = 43, 45, and 47° (circles, squares, and crosses, respectively). For large r n (and 
so for large AT) the line of sight crosses the open volume nearly centrally (/3 X <C 1.5^ v , 
T ~ 0), so that F ~ 0.1, regardless of the value of Cobs and a (see Fig. 3). Therefore, all the 
three numerical solutions stay close to the analytical solution for F — 0.1 (thin solid line). 
For smaller r n , (and AT < 1°), the numerical results diverge from each other: the case of 
Cobs = a = 45° (squares) remains close to the analytical solution with F ~ 0.1 (it follows 
eq. (19) with F & 1.07), because the parameter r oc /3 X is fixed and equal to zero. 

In the case of Cobs = 47° (equatorward viewing, crosses in Fig. 5), the parameter r increases 
with decreasing r n because the line of sight traverses more peripherally through the open 
volume. As can be inferred from Fig. 3 (panel for a = 45°) this makes F increase through 0.2 
up to ~ 0.26 for r ~ 1, and accordingly, the numerical solution in Fig. 5 crosses the dotted 
line for F = 0.2 and approaches the vicinity of the dashed line for F = 0.3. At logr n ~ —3.33 
the line of sight just grazes the outer boundary of the open volume (r ~ 1.08). At smaller 
radial distances the line of sight does not penetrate the open volume. 

In the poleward case of Cobs = 43° (circles in Fig. 5) , r becomes more negative with decreasing 
r n . Since the backward displacement of the open volume is stronger on its poleward side 
(see Figs. 2 and 3) the solution crosses the analytical curve for F = 0.2 (dotted line in 
Fig. 5) earlier than in the equatorward case (ie. at smaller |r|). At logr n ~ —2.43, the line 
of sight starts to cut the radiation beam above the notch visible in Fig. 2. This results in 
a discontinuous increase of F from ~ 0.25 up to ~ 0.39 (Fig. 3). Therefore, the numerical 
solution jumps to the vicinity of the dot-dashed line of F = 0.4 (Fig. 5). For more peripheral 
traverses (ie. for smaller r n and more negative r), F changes little between 0.39 and 0.35, 
and the numerical solution departs only slightly from the F = 0.4 curve. For logr n < —3.25 
the line of sight misses the open volume. This minimum value of logr n differs slightly 
from the one for /3 X = +2° (crosses) because in addition to the backward displacement, 
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the open volume is also slightly shifted (with respect to fl) towards the rotational equator 
(ie. downwards in Fig. 2). 


3.1.1. Determination of emission radius 

The above-described complicated form of the function AT(r n , a , Cobs) does not allow an 
easy determination of r n . Numerical determination of r n based on the known values of AT. 
a, and Cobs is rather sophisticated. Therefore, in Appendix B we discuss the particular cases 
when easy derivation of r n is possible and then we propose a procedure, which allows the 
possible range of r n to be constrained in a general case of arbitrary viewing geometry. Fig. 4 
and the examples considered in Appendix B make it clear that for a given phase shift AT, 
the misalignment formula predicts larger emission radius than the delay-radius relation of 
BCW91, and is therefore less prone to the problem of emission from the closed magnetic 
field lines. Another effect, which works in the same direction, ie. increases r n and may 
additionally lessen the problem is the altitude-dependent position angle swing. 


4. Altitude-dependent position angle swing 

When different parts of the pulse profile originate from different altitudes, the PA curve 
can no longer be described by the standard equation of Radhakrishnan Sz Cooke (1969). 
The simple analytical equation for the altitude-dependent PA swing is given in eq. (16) of 
BCW91. The formula predicts that if the radial distance r of the radio emission is uniform 
across the pulse profile, the entire PA curve is shifted rightwards (toward later phases) by 2 r 
with respect to the zero phase. This effect is illustrated in Fig. 6: The thin solid line with 
dots is the PA curve for the fixed emission altitude of r n = 10~ 2 , calculated numerically for 
the retarded dipole field with a = 45° and £ = 43°. In agreement with BCW91, its center 
lags the phase zero by 2 r n rad = 1.14°, and there is no discernible sign of the sweepback 
effect (see Section 2). 

If the central parts of the pulse profile originate from much lower radial distance r than 
the edge, and if r < 0.01R lc , the central parts of the PA curve do not exhibit any appreciable 
shift and nearly follow the undisturbed S-curve of Radhakrishnan & Cooke. The solid line 
which nearly passes through the (0, 0) point in Fig. 6b has been calculated for the emission 
from the last open magnetic field lines of the retarded dipole for the same a and £ as above. 
The corresponding radial distance of the radio emission as a function of phase <f> is shown in 
panel a of Fig. 6. Within the central parts of the PA curve, the emission altitude is negligibly 
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small and the PA curve nearly follows the undisturbed S-swing (dotted line in Fig. 6b). 

Thus, in the case when the central parts of the profile are emitted from much lower 
altitudes than the outer edge of the profile, the PA swing practically does not undergo the 
delay by 2r n . The misalignment formula (19) then becomes 

AT « 2 r n - Fr l J 2 (20) 

and is illustrated in Fig. 7 with the same layout as Fig. 4. One can see that in the considered 
case, the shift is much smaller than the shift predicted by the standard delay-radius formula 
(eq. 2, thick dashed line in Fig. 7). For a typical phase shift AT ~ 0.3° eq. (2) underestimates 
r by an order of magnitude. The range of altitudes for which the shift is negative is much 
larger than in the case with the constant emission altitude (eq. 19, Fig. 4). 


5. Conclusions and Discussion 

The rotational distortions of the vaccum dipole have twofold effect: in addition to the 
small changes of the local direction of the magnetic field, the region of the open magnetic 
field lines undergoes a strong distortion. 

The change of the local direction of B is second order in r n . We find, however, that 
it is symmetrical with respect to the (Q, ft) plane. The largest asymmetrical change of B 
direction is much smaller - third order in r n , in agreement with Shitov’s estimate. 

The region of the open field lines of rotating dipole possesses inherent asymmetry with 
respect to the (Q, ft) plane, at least in the limit of negligible plasma density. Therefore, the 
pulse window' is shifted tow r ard later phases with respect to the center of the position angle 
curve. The shift has the magnitude of the order of 0.2ry 2 . The open volume shift modifies the 
delay-radius relation if the center of the pulse profile is determined as the midpoint between 
the outer edges of the pulse profile, which are assumed to lie close to the outer boundary of 
the open volume. At low altitudes, where effects of aberration and propagation time delays 
are small, the open volume shift dominates and may result in the center of the PA curve 
preceding the center of the profile (negative AT). A majority of pulsars exhibits positive 
AT. The misalignment formula is consistent with this observation if the radio emission 
altitudes at the outer edge of their pulse profiles exceed ~ 10 -2 ili c (cf. Fig. 4) which is not 
an especially stringent constraint. 

The radii derived with the misalignment formula exceed those derived with the original 
delay- radius relation by a factor which increases quickly with decreasing altitude (and AT). 
This explains the trend of the delay-radius relation to predict emission radii smaller than the 
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geometrical radii. The underestimate may also be produced/enhanced by the low-altitude 
emission within the central parts of the pulse profile. When both these effects work together, 
the standard delay-radius relation may underestimate r by an order of magnitude even for 
the relatively large phase shifts (A$ ~ 0.3°). 

The influence of the open volume shift on the method based on the core-cone shift 
(eq. 3) is difficult to assess, because the locations of the conal maxima do not need to follow 
the outer boundary of the open volume. If they did, eq. (3) would have to be replaced with 
A<I>drh ~ 2r n - Frl/ 2 . If the cones have the shape of the open boundary, but their size is 
smaller, the backward shift is decreased and equals A^ ov ~ F- m (9- m /6 b )r D , where 9- m < 9 h is 
the opening half-angle of the cones, and Fj n = F(t = T; n ) with T- m — fi x /9- m . Eq. (3) in this 
case becomes A<f> DRH « 2 r n - F m {9- m /9 h )r l J 2 . The influence of the open volume shift on the 
geometrical method is small because the rotation increases the transverse size of the open 
volume insignificantly (by a factor smaller than ~ 1.2, see Fig. 2). 

Our misalignment formula suffers from the same internal inconsistency as the theory of 
BCW91: On one hand relativistic electrons are assumed to participate in the corotational 
motion of the magnetosphere, which therefore must be filled in with the Goldreich- Julian 
charge density. On the other hand the vacuum structure of the magnetic field is assumed 
(BCW91 used the vacuum magnetic dipole of static shape, see their eq. 9). Like the sweep- 
back effect, a current flow may affect the delay-radius relation (as well as our misalignment 
formula) in two ways: either by the local distortions of the magnetic field at the emission 
point, or by the displacement of the open volume boundary as determined by the conditions 
near the light cylinder. Locally, the currents are a higher order effect than considered in this 
paper and we neglect them (see Hibschman & Arons 2001). 

As we have seen on the example of the sweepback, this does not need to necessarily imply 
that the currents near the R \ c have negligible impact on the outer boundary of the open 
volume. Since there is no well established theory of plasma loss near the light cylinder, we 
are forced to neglect its influence on the geometry of the open volume (as did BCW91). 
With two incomplete theories at hand (ie. BCW91 and our), a reader may be more inclined 
to use the delay-radius relation of BCW91, as the simpler one. However, as we have shown 
in Section 2.1, the assumption of BCW91 about the symmetry of the open volume is incor- 
rect at least in the limit of the negligible plasma density. Whereas we cannot prove that 
the plasma effects nea,v-R\ c are negligible, it is unlikely that they completely cancel out the 
asymmetry generated by the rotation. 

We thank U. Dyks for the derivation of eq. (9). JD thanks B. Rudak for all the years of 
fruitful collaboration. This work was performed while JD held a National Research Council 
Research Associateship Award at NASA/GSFC. This work was also supported by the grant 
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PBZ-KBN-054/P03/2001 (JD) and by the NASA Astrophysics Theory Program (AH). 

A. Retarded vs static dipole 

In the reference frame with z || Q, the “retarded” field B TBt of the vacuum magnetic point 
dipole with the magnetic moment given by eq. (5) has the following cartesian components: 

-Bret,* = ^ 1 3xz cos a 4- sin a ( [(3x 2 - r 2 ) + 3 xyr n + (r 2 — x 2 )r 2 ] cos (fit — r n ) + 


[3xy - (3x 2 - r 2 )r n - xyr 2 ] sin(f2t - r n ) ) } (Al) 

, y = ^ {3yzcosa + sina ( [3 xy + (3y 2 — r 2 )r n - xyr 2 ] cos(Qt - r n ) + 

[(3y 2 - r 2 ) - 3 xyr n + (r 2 - y 2 )r 2 ] sin(ftt - r n ) ) } (A2) 

B ret , 2 : = ^ { (3z 2 — r 2 ) cos a + sin a ^ [3xz + 3 yzr n — xzr 2 ] cos (fit — r n ) 4- 

[3 yz - 3 xzr n - yzr 2 ] sin(ftt - r n ) ) } (A3) 

It can be derived by any of the methods described in Yadigaroglu (1997), Arendt & Eilek 
(1998), or CRZ2000. Also, it can be obtained by taking the limit of Ra S /R\ c C 1 in the 

solution of Deutsch (1955). For any position r 0 , time to and time interval At it holds that 

Bret (to 4- At, R z (OAt)f 0 ) = R z (QAt)Rret (to, r 0 ) , (A4) 


where R z (u)f represents the rotation of the vector r by the angle u around the z axis. Thus, 
the dependence on the time t only reflects the rigid rotation of the pattern of distorted 
magnetic field lines around the z axis with the angular velocity 0 z. The magnetic field B st 
of the static-shape dipole associated with the same magnetic moment (2 (eq. 5) is given by 
the same set of equations (Al - A3) with r„ = 0. 

Using eqs. (Al - A3) with t = 0 it can be immediately shown that the difference 
A B — Bret — Bst is given by eqs. (6 - 8). Using (4) and (6 - 8) one can find that the angle 
k — arccos(B ret • B st / (B ret Bst)) is given by eq. (9) with f u / 2 , and f 3 given by 

fi = h- y? , f 2 = h 2 ~ 2g x g 2 (A5) 

(A6) 

(A7) 


and 

with 


fz — 2 gf +h 3 — g 2 — 2g x (g 3 4- h x ) 
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where J5 st)X is given in eq. (13). 


The cartesian components of the retarded magnetic field given by eqs. (A1 - 
rewritten into the following spherical components of B Tet in the reference frame 


(A8) 

(A9) 

(A10) 

(All) 

(A12) 

(A13) 

(A14) 


A3) can be 
with z II Q: 


B r et ,r = 

2n 

— {cos a cos 9 + sin a sin 9 [r n sin A + cos A] } 

(A15) 

■Bret ,</> ~ 

— -^-sina [(r„ — l) sin A + r n cosA] , 

(A16) 

Bret,0 = 

^ {cos a sin 9 + sin a cos 9 [— r n sin A + (r\ — l) cos Aj } 

(A17) 


A = r n + <j) — fi 1. 

(A1S) 


where 

Again, the (spherical) components of B st are given by eqs. (A15 - A18) with r n = 0. 
Using eqs. (A15 - A17) in eq. (4) it can be easily shown that 
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r l sin S - ^ rl cos 6 + O (r*) 
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(A19) 

(A20) 

(A21) 


where 5 — <f> — Clt. 
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B. Determination of emission radius 


As can be inferred from Fig. 3. the value of F is (nearly) fixed and close to 0.1 whenever 
|r| <C 0.4, ie. whenever \f3\ 0.6ry 2 . The condition is fulfilled for any r n and a when fi « 0 

(a ~ Cobs, squares in Fig. 5). The other case when the condition is fulfilled is when the 
measured shift between the position angle center and the profile center is large (> 1°) and 
the impact angle \/3\ is small: the radiation comes then from high-altitudes, where » m. 
This is why all numerical results shown in Fig. 5 approach the solid line of F = 0.1 when 
A<f> increases above 1°. 


When F is fixed (like in the above-described cases), the equation (19) can be inverted 
to obtain analytical solutions for r n (A$): 


r n = 32" 1 


8A$ + F 2 + k (16A $F 2 + F 4 ) 1/2 


(Bl) 


where the parameter k = 1 for A$ > 0 whereas k — ±1 for A$ < 0. Thus, the solution for 
r n is unique if A$ is positive, whereas for the negative A$ two solutions are possible. One 
can attempt to reject the smaller one of these solutions by using the theoretical constraints: 

R 4 

r n > and r n > - ft 1 . (B2) 

The second constraint holds when the emission is limited to the open volume. 3 In the two 
above-mentioned cases (when « 0 or A<£ > 1°), one can directly use eq. (Bl) with F = 0.1 
to calculate r n . 


In a general case (including the common case when f} ^ 0 and the shift is small) one can 
constrain r n to a degree w 7 hich depends on the information available. Without the knowledge 
of a and /3 one can obtain a rough estimate of r n based on the measured shift A$ by assuming 
some averaged, fixed value of F in eq. 19 (eg. F = 0.2). This is equivalent to igoring the 
dependence of F(a, r) shown in Fig. 3. For r n < 10 -2 this may give results wrong by as 
much as an order of magnitude due to the large horizontal spread of r a as allowed by the 
range of F (Fig. 4). 

Therefore, a better approach is to use the following two-step procedure: 1) Using eq. (Bl) 
one can calculate the range of r n allowed by 0.1 < F < 0.4 and by the conditions (B2). Or, 
one can draw a horizontal line in Fig. 4 to determine graphically the range of r n allowed by 
the range of F (the measurement error of A$ can easily be taken into account by drawing 


3 In application to the retarded dipole field this condition is approximate, not exact. Therefore, the 

minimum values of log r n determined numerically for $ = ±2° in Fig. 5 (circles and crosses) differ slightly 

from —3.266. 
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a horizontal strip instead of the line). For —0.0358° < AT < 0, two ranges of r n are 
allowed and again the constraints (B2) may be used to try to reject the lower range. For 
—0.573° < AT < -0.0358°, the allowed range of r n is between r™ in and r™ ax where r™ ,n is 
given by eq. (Bl) with k = — 1 and F — 0.4 and r“ ax by eq. (Bl) with k — +1 and the same 
value of F = 0.4. 

2) With the knowledge of a and j3 one can more tightly constrain the allowed range of r n 
by calculating r m ; n = r(r® ax ) and r max = r(r™ in ) where r“ in and r™ ax were determined in 
the previous step. Then, from Fig. 3, using the appropriate panel, one can determine the 
narrower range of (F m i n , F max ) corresponding to the range of (r m ; n , r max ). Then, one returns 
to the step 2) in which the tightened range (F m i n , F max ) must be used instead of the original 
range of (0.1, 0.4). 

If the sign of fi is unknown, the step 2) may (but does not have to) give two ranges of 
r n . If the sign is positive (equatorward viewing) the range of F considered in the step 1) 
may be narrowed to 0.1 — 0.25 (Fig. 3). 

As an example we apply the method to PSR B0301+19 and B0525+21. Although the 
pulsars exhibit very simple pulse profiles and smooth PA swings, they are among the most 
problematic objects for the delay-radius relation, which places their radio emission well into 
the region of the closed magnetic field lines (see Table 3 in Dyks et al. 2004a). In the case 
of B0301+19 ( P = 1.38 s), BCW91 find AT = 0.2 ± 0.1° at 1.4 GHz. The range of r n 
for F between 0.1 and 0.4 is (r“ m ,r“ ax ) = (1.4 • 10 -3 ,1.25 • 10 -2 ) with the error of AT 
included (step 1). For /3 X = —0.96 ± 0.63° (Everett & Weisberg 2001, hereafter EW2001) 
Fig. 3 gives (r min , r max ) = (—0.5, —0.03) which includes also the error of /3 X . EW2001 derived 
a — 162.4 ± 11.8°, so that tx — a « 17.6° and we may use the panel of Fig. 3 for a — 20° 
to constrain the range of F to (F min , F max ) = (0.14,0.22). This new range of F translates 
to (r™ m ,r“ ax ) = (2.0 • 10“ 3 , 5.3 • 10~ 3 ) (with the errors of both AT and /3 X included). None 
of the constraints (B2) narrows this range. This result is in agreement with the condition 
r n > r geo /F lc . For the observed pulse width W « 15.9° (BCW91) and for a and j3 cited 
above, one finds r geo /i?i c ss 1.3 • 10~ 3 . Our value of r n = (2 — 5.3) • 10“ 3 implies that 
the emission associated with the profile edge must originate from magnetic field lines with 
magnetic colatitude 6 m « (0.5 — 0.8)#£ v . The original delay-radius relation gives r n 9- 10~ 4 
(3 times smaller) which results in 9 m fd T oy > 1 (emission from the region of closed field lines). 

In the case of B0525+21 (P = 3.74 s), BCW91 find A$ = 0.3±0.1° at 430 MHz. At this 
value of A$ the range of F = 0.1 — 0.4 translates into (r™ m , r™ ax ) = (2 • 10~ 3 , 1.3 • 10 -2 ) with 
the error of A<3> included. For /3 X = 1.5 ±0.08° (EW2001) we find (r m j n ,r max ) = (0.15,0.39). 
EW2001 suggest i r — a = 63.2° which implies that F ~ 0.1 (Fig. 3, panel for a = 60°). 
Using this value in eq. (Bl) one finds that r n = (2.6 ± 0.6) • 10" 3 which is two times larger 
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than the value predicted by the original delay-radius relation (eq. 2). However, the value 
is still six times smaller than r geo /R\ c = 1.7 • 10“ 2 (calculated for W = 20.4°) and implies 
that the emission at the profile edge comes from magne tic field lines w a o j. a tudes 0 m — 
(2.4 — 3.O)0£ V , ie. from the closed line region. A large part of this discrepancy can be removed 
by considering different emission altitudes across the pulse profile. 

In the case of the low-altitude emission within the profile center (considered in Section 
4), the determination of the radial distance of radio emission r n is performed as above, with 
the only difference in that eq. (Bl) should be replaced by the inverse of eq. (20), ie. Fig. 7 
should be used instead of Fig. 4. Even allowing for the low-altitude origin of the central 
parts of the pulse profile, for B0525+21 at 430 MHz we still have r n < r geo /R] C . Apparently, 
either the radio beam of this pulsar has irregular shape, or other effects are important. 
These may include the overestimate of dipole inclination and/or |/?|, the broadening of the 
observed pulse width due to the low energy of radio emitting electrons (all of which lead to 
the overestimate of r geo ) or refraction, eg. Lyubarski & Petrova (1998); Fussell & Luo (2004). 
However, our radial distance r n ~ 7- 10~ 3 underestimates r geo /Ri c by a factor of 2.5 whereas 
the original delay-radius formula underestimates it by a much larger factor of 13. 
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Fig. 1. — The influence of rotation on the local directi on of the magnetic field. The vectors 

of B ret and B st are shown for two points in the equatorial plane of the orthogonal rotator 

(a = 90°). The points are located symmetrically with respect to the (fi, p) plane (the plane 

is orthogonal to the page and contains x). The rotation is to the right and it is assumed 

that x » y (region near the dipole axis) and x <C R\ c - With accuracy of r„, the retarded 
— * 

field B Tet is symmetrical with respect to the (Q,p) plane, ie. /q = /q and \fa\ — <j> \. More 
precisely, /q = K t + Q{rl) and <f> t = -fa - O(r^), ie. K t < and \fa\ > fa. 
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Fig. 2. — Crossections of the open volume with the star-centered sphere of radius r <C R\ c 
calculated for the retarded dipole field B Tet (thick solid). The circles have radius of rr„ /2 and 
are centered on the magnetic moment /2 which protrudes perpendicularly from the page at 
the center of each panel [(a; m , y m ) = (0, 0)]. The magnetic field B ret which permeates each of 
the panels is nearly the same as that of the static-shape dipole with the axis containing p. 
Rotation is to the left. The backward displacement of the retarded contours with respect to 
the (0, 0) points (or circles) results in the shift of the center of the pulse profile toward later 
phases. The shape of the contours does not depend on r and P as long as r < Ric ( P = 1 
s, and r n = 0.01 was assumed in the figure). Their size scales as rrU 2 . The dotted ovals 
present the case of the static-shape dipole. 
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Fig. 3. — The parameter F of eq. 17 as a function of the parameter r = /? x /(l .hr]/ 2 ) for the 
same dipole inclinations a as in Fig. 2. Though the results were obtained for r n = 0.01, the 
curves change little with r n , as long as r n < 1. 
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Fig. 4. — The phase shift between the center of the PA curve and the profile center as a 
function of the radial distance of the radio emission. The four curves show the case with the 
rotational distortion of the open volume included (eq. 19) and correspond to different values 
of the parameter F shown in the lower right corner of the upper panel. The thick dashed 
line presents the original delay-radius relation (eq. 2) which does not include the sweepback 
effect. In the upper panel the positive shift is shown (the PA curve lags the profile). The 
lower panel is for the negative shift (the PA curve precedes the pulse profile). Note that the 
original delay-radius relation significantly underestimates r n , especially for small and 
r/Ric- 




Fig. 5. — Three numerical results for a = 45° (and different viewing angles £) overplotted 
on the curves from Fig. 4. The circles are for the poleward viewing geometry with £ = 43° 
(the curve is broken into two parts, one with positive, and the other with negative A<$), the 
crosses are for the equatorward viewing with £ = 47°, and the squares are for a = £ = 45°. 
For more details see text. 
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Fig. 6. — Influence of variations of emission altitudes across the pulse profile on the shape 
of the PA curve. Thin solid line with dots (panel b) presents the position angle curve for 
emission from the fixed radial distance of r n = 0.01. Its center lags the phase zero by 
2 r n rad « 1.14°. The thick solid line in panel b is for the radio emission from the last open 
field lines, ie. it corresponds to different radial distances shown in panel a. Note that the PA 
curve for the case of varying r does not exhibit any noticeable lag. The spread in r visible 
in panel a) near 4> = — 7° corresponds to the notch in the open volume which appears for 
moderate dipole inclinations (see Fig. 2, a = 45°). The dotted line (nearly overlapping with 
the thick solid) is the curve of Radhakrishnan & Cooke (1969), undisturbed by the special 
relativistic effects. 
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Fig. 7. — Dependence of the shift between the center of the PA curve and the profile center 
on the radial distance of the radio emission in the case when the central parts of the pulse 
profile originate from much lower altitude than its edge (ie. for the case marked in Fig. 6 
with the thick solid line). The layout is the same as in Fig. 4. Note the increased divergence 
from the original delay-radius relation (thick dashed line). 




